Superconductivity of MoaSby from first principles 
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Superconductivity in MoaSby is analyzed using the combined electronic structure and phonon cal- 
culations, and the electron-phonon coupling constant Aph = 0.54 is determined from first principles. 
This value explains the experimental value of the superconducting critical temperature Tc — 2.2 K. 
The possible influence of spin fluctuations and spin gap on the superconductivity in MoaSbr is 
discussed, and electron-paramagnon interaction is found to be weak. 
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A paramagnetic intermetallic compound MoaSbr is 
a type II superconductor,^-^ with the critical temper- 
ature Tc ~ 2.2 K. The temperature characteristics of 
the specific heat, the superconducting gap, and the 
magnetic critical field suggest that the conventional 
electron-phonon interaction might be responsible for the 
superconductivity.— '^i^'^ Recently, however, Candolfi et 
al.^ argued that spin fluctuations (SFs) are present in 
MosSby. This interpretation is supported by two un- 
usual features: (i) the quadratic temperature dependence 
of both electrical resistivity and magnetic susceptibility, 
as well as (ii) the high value of the susceptibility at room 
temperature. They also reported a much smaller value of 
the electronic specific heat jump^ at the transition point 
AC/^Tc = 1.04 than the weak-couphng BCS value 1.43, 
which might suggest additional enhancement of the elec- 
tronic specific heat coefficient by the SFs. Very recently, 
Tran et al^ observed a peak in the specific heat Cp{T) 
at T* = 50 K, which was interpreted as supporting the 
presence of spin gap. Also, they explained the anomalous 
behavior of the magnetization and resistivity in terms of 
the gap opening. Moreover, they analyzed the electronic 
specific heat in the superconducting state in terms of 
the two BCS gap model,'* and reported a higher value of 
AC/'jTc — 1.56 than the one measured beforei^ 

In order to elucidate the possible origin of supercon- 
ductivity, an ab initio approach which involves the elec- 
tronic structure and phonon calculations may be used to 
determine the electron-phonon coupling (EPC) constant. 
For instance, a recent determination of the EPC con- 
stant suggested that the superconductivity in PuCoGas 
is driven by an unconventional mechanism based on an- 
tiferromagnetic (AF) fluctuations.^ Here we present an 
ab initio study of the EPC constant and superconduc- 
tivity in Mo3Sb7, where SFs might play a role. The 
electron-phonon interaction is treated within the rigid 
muffin tin (MT) approximation. The superconducting 
critical temperature Tc and its possible modification by 
SFs is discussed using two approaches: (i) the McMillan 
formulaf^i^ and (ii) the equation for Tc which explicitly 
includes the presence of paramagnons.^^ 

Electronic structure calculations were performed using 
the Korringa-Kohn-Rostoker (KKR) multiple scattering 



method. "'^"'^ The crystal potential was constructed in the 
framework of the local density approximation (LDA), us- 
ing von Barth and Hedin formulcfci^ for the exchange- 
correlation part. For all atoms angular momentum cut- 
off Imax = 4 was set; k-point mesh in the irreducible part 
of the Brillouin zone (BZ) contained about 400 points. 
Density of states (DOS) was computed using the tetra- 
hedron k-space integration technique, generating about 
1500 tetrahedrons in the irreducible part of the BZ. 
Semirelativistic calculations results are presented here. 
Since our main goal in this work is to estimate the EPC 
constant from first principles within the rigid MT approx- 
imation, spherical potential approximation for the crystal 
potential is used, as is required in this approach. MosSby 
crystallizes in a cubic bcc structure (space group ImSm) 
of the IrsGey type, with lattice constanti2, a — 9.58 A. 
The primitive cell of MoaSby contains two formula units, 
i.e. 20 atoms, occupying three nonequivalent positions: 
Mo in (12e) with x = 0.3432, Sb(l) in (12d) and Sb(2) 
in (16f) with x = 0.1624. 

The phonon frequencies were determined within the di- 
rect method^ which utilizes Hellmann-Feynman forces 
obtained by performing small atomic displacements of 
nonequivalent atoms from their equilibrium positions. 
From them the dynamical matrix is determined and 
diagonalized to obtain the phonon frequencies at each 
wave vector. The crystal structure optimization and 
calculations of the complete set of Hellmann-Feynmann 
forces were performed using the first-principles VASP 
package^^ which makes use of the Perdew, Burke, and 
Ernzerhof (PBE) functionalii^ The calculations were per- 
formed on a X x 1 supercell (containing 80 atoms) 
with periodic boundary conditions. The wave func- 
tions were sampled according to Monkhorst-Pack scheme 
with a k-point mesh of (4,4,4). After the optimization 
we obtained the lattice parameter a = 9.6405 A and 
the atomic positions of (0.3421,0,0), (0.25,0,0.5) and 
(0.1608,0.1608,0.1608) for Mo, Sb(l) and Sb(2), respec- 
tively. The determined values are in very good agreement 
with the experimental dataj^ 

The electronic structure was computed for the exper- 
imental lattice parameters and atomic positions. Total 
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and site-decomposed electronic DOSs of MosSby are pre- 
sented in Fig. [TJ The most intriguing feature of the elec- 
tronic spectrum is the presence of a narrow band gap 
just above the Fermi level, with Ep located in the range 
of sharply decreasing DOS.^^ In the inset in Fig. [T] one 
observes that Ep coincides with a local DOS maximum. 
By analyzing the angular contributions to the total DOS 
at Ep, presented in Table [H we deduced that the bands 
near Ep are built out of the Mo{4d) and Sb(5p) states. 
The largest atomic contribution comes from Mo atom, 
with the value nuoiEp) ~ 14 Ry~-'^/spin, being not far 
but below the magnetic instability [the computed Stoner 
parameter / satisfies InMo{Ep) ~ 0.7]. Note that the 
spin-polarized KKR calculations assuming ferromagnetic 
(FM) spin order led to the nonmagnetic ground state. 

Tran et al.^ suggested the opening of the spin gap be- 
low 50 K, caused by the AF interactions between the 
selected nearest pairs of Mo atoms. They argued that 
these atoms form dimers, and the AF interaction sta- 
bilizes there spin singlets (but long-range order is ab- 
sent). We examined a few possible AF structures for this 
compound, e.g. with alternating moments in Mo planes, 
but stable AF configuration could not be reached and 
all magnetic moments converged to zero values. Note, 
that the proposed model;^ including one AF and two FM 
types of Mo-Mo interactions, creates a geometrical frus- 
tration of the Mo sublattice. The high value of the DOS 
at Ep, as well as the suggested different magnetic inter- 
actions between Mo atoms, may also give rise to the SFs, 
which could appear in real sample. 

The electronic structure results were used to cal- 
culate the electronic part of the EPC constant, i.e. 
the McMillan-Hopfield r]i parameter a^i'^^ for each atom. 
They follow from the formulai ^^'^" 
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FIG. 1: (Color online) Total and site-decomposed densities 
of electronic states in MosSby (per formula unit). The inset 
shows the details of the DOS near Ep = 0. 



where V{r) is the self-consistent potential at site i, i?MT 
is the radius of the i-ih. MT sphere, Ri{r) is a regular 
solution of the radial Schrodinger equation (normalized 
to unity inside the MT sphere), ni{Ep) is the ^~th par- 
tial DOS per spin at the Fermi level Ep, and N{Ep) is 
the total DOS per cell and spin. The values of rji pa- 
rameters ([l}, with contributions from each ^ — > Z -f 1 
scattering channels, are presented in Table HI For Mo, 
the d-f channel is the most important one (typically for 
d-element), whereas p-d contribution dominates for both 
Sb atoms. The Sb(l) and Sb(2) atoms have very similar 
rji parameters, despite quite different p-DOSs. This is a 
result of opposite behavior in both partial DOSs, i.e. for 
the Sb(2) atom the lower p-DOS is compensated by the 
larger value of d-DOS [the radial wave functions matrix 
elements form Eq. ([T|) are similar in both cases]. 

The phonon dispersion relations along the high symme- 
try directions and the total and site-decomposed partial 
phonon DOSs were computed by random sampling of the 
BZ and are presented in Fig. [5] for the optimized super- 
cell. The optic phonons give three characteristic maxima 
of the phonon DOS p{u)) at a; ~ 2.8, 4.4, and 6.3 THz. 
The Mo atoms, which are about 30% lighter than Sb 
atoms, contribute mainly to the high frequency part of 
the phonon DOS. The phonon DOS was used to compute 
the average square site-decomposed phonon frequencies 
{bjf) presented as well in Table [D These quantities, to- 
gether with {rii\ parameters, are needed to deduce the 
EPC constant 

I '-LI J 

Here i runs over all the atoms in the primitive cell and Mi 
is the atomic mass. For a review, more detailed discus- 
sion of the approximations involved in this approach, and 
a number of references to the previous rigid MT studies, 
see e.g. Ref. j21) and references therein. 




FIG. 2: (Color online) Phonon dispersions along the high 
symmetry directions of the BZ (left) and total and site- 
decomposed densities of phonon states in MoaSby (right). 
The special points are: T = (0,0,0), H = (1/2,-1/2,1/2), 
P = (3/4, -1/4, 3/4), r = (1, 0, 1), N = (1, 0, 1/2). 



3 



TABLE I: Site-decomposed electronic and dynamic properties of MoaSby. ni{EF) is in Ry ^/spin, rji in mRy/as^ (both per 
atom), u}i in THz. Values of \i take into account the number of i-type atoms in the primitive cell: 6 Mo, 6 Sb(l), 8 Sb(2). 
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Surprisingly, one finds that all the atoms are equally 
important for the onset of superconductivity in Mo3Sb7. 
The contribution from Mo atoms to the total Aph, despite 
the dominant character of Mo states near Ep, is only 
slightly larger than those from Sb(l) and Sb(2) respec- 
tively. This is a consequence of higher partial phonon 
frequencies for Mo. It is worth noting that Sb(l) and 
Sb(2) have the same Xi values in spite of rather different 
average phonon frequencies. Here, the effect of higher 
(uif) for Sb(2) is compensated by the larger multiplicity 
of this crystallographic site. The calculated total EPC 
constant ^ is Aph = 0.54, which quahfies MoaSbr as the 
medium-coupling superconductor. 

We estimated the superconducting critical temper- 
ature Tc using two formulas: (i) a McMillan-type 
formula^i^ with the logarithmically averaged phonon fre- 
quency Wph = (wiog) in the prefactor. 



T, 



Wph 

r2o 



exp 



1.04(1 + Aeff) 



Aeff-Aieff*(l + 0.62Aeff) 



(3) 



and (ii) the formula including the interaction of elec- 
trons with paramagnons, and successfully applied before 
to MgCNi3,io 



Tr = 1.14 w 
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, (4) 



(5) 



Here Asf stands for the electron-paramagnon interaction 
parameter, and Ugf is the characteristic SF frequency 
(temperature) . 

The interplay between SFs and superconductivity is a 
well-known problem in the theory of superconductivity. 
In conventional superconductors, with electron-phonon 
pairing mechanism, FM SFs are known to compete with 
superconductivity, leading e.g. to the lack of super- 
conductivity in palladiumi^ More recently, SFs (para- 
magnons) were studied in the context of superconductiv- 
ity in MgCNi3,^'^ or for elemental metals under pressure: 
Pg^ 23,24 g^j^jj Sci^ In fact, one finds that in case of SF su- 
perconductor the McMillan formula^'^ may still be used, 
but the parameters A and fi* , applied when Agf = in 
Eq. ([31), are then renormalized tOf2? Agff — Aph/(1 + Asf), 

Meff* = (m* + Asf)/(1 + Asf). 



First we calculate Tc without taking into account the 
SFs, i.e. using Eq. ^ with Aeff = Aph and /^eff* = M*- 
Since the value of Coulomb pseudopotential parameter 
fi* is unknown, we present Tc in a realistic range of 
0.08 < ^i* < 0.18 in Fig. H For the typical values of 
II* and the calculated Wph = 143 K we get Tc = 2.4 K 
ill* = 0.10) and 1.6 K (/.t* = 0.13). Note, that when the 
prefactor in the McMillan equation is set to the origi- 
nal value^ 0/1.45, and the experimental value'^ of Debye 
temperature = 310 K is used, the resulting temper- 
atures are higher: Tc = 4.0 K {fi* = 0.10), 2.6 K {^i* 
— 0.13), 1.8 K {fi* = 0.15). These results demonstrate 
that, depending on the prefactor, the experimental criti- 
cal temperature Tc = 2.2 K may be explained using the 
EPC constant Aph = 0.54 derived within the rigid MT 
approximation, and taking fi* between 0.10 and 0.13. 

Next we analyze the possible influence of SFs on the 
transition temperature Tc. The electron-paramagnon 
mass enhancement Asf is treated as a parameter. It is 
important to note that if one explicitly takes into account 
the SF effect on the superconductivity, the starting value 
of fi* (i.e. before its renormalization by Agf ) can be taken 
smaller tha n ty pically used (e.g. for Nb /x* — 0.086 was 
used in Ref.|2y). Since Eq. ([4]) involves additional param- 
eter, i.e. the characteristic paramagnon frequency Wgf, in 
this case we plotted Tc against LOgf for some representative 
values of Asf in Fig. [H For Wgf > 100 K the value of Tc is 
practically nonsensitive to the chosen ajgf, thus this value 
was used in the calculations. Using Eq. one finds that 
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FIG. 3: (Color online) Critical temperature Tc as a function 
of Coulomb parameter for Asf — and two prefactors in 
Eq. ©: Wph/1.20, and 9/1.45. Parameters: Aph = 0.54, 
Asf = 0, Wph = 143 K, e = 310 K. 



4 



3 



i 1 i 1 1 i 
^.^ 

\ 


1 1 i 1 i 1 1 i i 

Asf=0.03 

Asf=0.05 


Asf=0.08 " 

Asf=0.10 ■ 








, . , 1 . , 


. 1 , . , 1 . , , 





50 100 150 200 250 300 

a;,f(K) 

FIG. 4: (Color online) Critical temperature Tc as obtained 
from Eq. Q for increasing paramagnon frequency ujst and for 
different values of Agf. Parameters: Aph = 0.54, = 0.09, 
cjph = 143 K. 

temperatures close to the observed Tc — 2.2 K may be ob- 
tained for Asf = 0.03 and fi* = 0.08 — 0.09 (corresponding 
to the effective fi*ff = 0.11 - 0.12), i.e. Tc = 2.3 K and 
Tc = 2.0 K, respectively. Fig. [4] shows that Tc quickly 
tends below 2 K when the electron-paramagnon inter- 
action parameter Asf > 0.05i^ Thus we conclude that 
the observed magnitude of the superconducting critical 
temperature can be explained taking into account the SF 
effects, but the Agf parameter has to be relatively small, 
Asf — 0.03, if the EPC parameter Aph = 0.54 obtained in 
our study is used. 

Another interesting question concerns the influence of 



the spin gap, detected below T* = 50 K, on the supercon- 
ducting state of MosSby. In view of the present results, 
this effect cannot be very strong, since (i) not all the Mo 
atoms are involved in building the singlet dimers (respon- 
sible for the gap^), and (ii) the Mo sublattice contribu- 
tion to the total EPC constant Aph is about 35%, with 
the rest provided by the two Sb sublattices. 

In summary, the results of electronic structure and 
phonon calculations were used to calculate the Aph pa- 
rameter for the spin-fluctuation/spin-gap superconduc- 
tor MoaSby, within the rigid MT approximation. The 
estimated value of Aph = 0.54 quahfies MosSby as a 
medium-coupling superconductor. The experimentally 
observed critical temperature Tc ~ 2.2 K may be cor- 
rectly reproduced even including the presence of param- 
agnons, with small Asf — 0.03. Thus, the spin fluctua- 
tions may exist in MoaSby, but the electron-paramagnon 
interaction has to be moderate. Since the Mo contri- 
bution to the constant Aph is estimated to be compara- 
ble to Sb(l) and Sb(2) sublattices, the possible influence 
of spin gap on the superconductivity is expected to be 
rather weak. However, in the range of the EPC constant 
Aph ~ 0.5 the value of Tc is quite sensitive even to small 
changes in Aph, so a more quantitative explanation of the 
superconductivity in MoaSby requires further study. 
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